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Cosmic ray protons propagating in a spatially-homogeneous but time-dependent field of axions 
or axion-like particles (ALPs) emit photons in a way that is reminiscent of Cherenkov radiation by 
charged particles in a preferred background. We compute the emission rate and energy spectrum of 
the photons, and discuss the possibility of their detection using the Square Kilometre Array which 
is currently under construction. In the case of a non-detection, constraints can be placed on the 
parameter space of ALPs whose mass lie between 10~^eV and 10~^ eV under the assumption that 
they are the primary constituent of dark matter. 

I. INTRODUCTION 

We are now in an exciting era for particle physics. In the high energy frontier, the Large Hadronic Collider has begun 
operation, giving us unprecedented access to physics at the TeV scale. What has perhaps garnered less attention, 
but is certainly no less important, is that time may also be ripe for the discovery of new physics at the sub-eV scale. 
Currently, a large number of experiments around the world are actively searching for the axion, its cousins axion-like 
particles (ALPs), and various light particles (e.g. chameleon, hidden photons, mini-charged particles, etc.), placing 
ever improving limits on their parameter spaces. 

The exploration of particle physics at the sub-eV scale was initially motivated more than thirty years ago by the 
proposed solution by Peccei and Quinn to the strong CP Problem, namely the absence of CP violation in the strong 
interaction [1] [5] . This proposal, which predicts the existence of a new pseudoscalar particle (the axion) , remains the 
most convincing solution to the strong CP Problem as of today. Since then, two developments have injected further 
impetus into the field of low-energy particle physics. First, it has been realized that the axion has exactly the right 
properties to be dark matter [3HS] • Secondly, spin-0 particles akin to the axion (known as ALPs) are found generically 
in string theory [7]. It is thus no wonder that axions and ALPs have become something that is being sought avidly 
by experimentalists all over the world. 

In the first axion model by Weinberg and Wilczek, the symmetry breaking occurs at the weak scale, though this was 
quickly ruled out by experiments [SI |H] . It was subsequently realized that the symmetry breaking scale could be much 
larger, though this implies that the axion couples extremely weakly to ordinary matter |10H13j . They have thus come 
to be known as "invisible" . Nonetheless, it has been demonstrated that their existence can be probed by exploiting 
their coupling to photons |14H16| . This idea has been implemented in various experiments around the world. The 
ADMX experiment is a realization of the concept of the axion haloscope, in which some halo axions in a magnetic 
field are induced to convert to microwave photons, which can then be picked up by an antenna jTT]- The CERN 
Axion Solar Telescope [18] and the Tokyo Helioscope [19j are axion helioscopes which convert axions from the Sun 
into X-rays in a magnetic field. On the other hand, photon regeneration experiments ( "shining-light-through- wall" ) 
do not rely on external sources of axions [20 . Instead, (some) photons in a laser beam propagating towards an opaque 
wall are converted into axions in a magnetic field. The axions travel unimpeded through the wall, behind which is an 
identical setup of magnets, where some axions are converted back to photons which can be detected. 

In this paper, we propose a new observational probe of axions or ALPs, based on the fact that charged particles 
propagating in a time-dependent axion field emit photons. There are three processes that contribute to the photon 
emission, which are shown in figures[l][2] and[3j In a sense, this phenomenon shares some resemblance with Cherenkov 
radiation that arises from spontaneous Lorentz violation, as the axion field effectively picks out a preferred rest frame 
in the universe (the one in which it is homogeneous). The time-dependent axion field can be viewed as a source 
of energy (due to its time dependence), but not three-momentum (due to its homogeneity). Consequently, some 
processes that are kinematically forbidden (energy and momentum conservation cannot be satisfied simultaneously) 
before now become possible. 

Understandably, because the axion's coupling to ordinary matter is very small, the rate of photon emission is 
extremely tiny. In fact, with current technologies, performing a laboratory experiment to observe the emission is 
definitely not feasible. A rough estimate shows that an electron accelerator with a reasonable length and electron flux 
will have to be in operation for more than the age of the universe before the emission of a single photon. Fortunately, 
this phenomenon occurs naturally in cosmology: by cosmic rays (primarily protons) propagating in a time-dependent 
axion field. As we will demonstrate, the abundance of axions (or ALPs) and cosmic rays in our galaxy compensate for 
the smallness of the coupling. With the aid of a detector with a collecting area of 10"'^'' cm^ (for example, the Square 
Kilometre Array currently under construction), this can give rise to a weak but detectable signal. 
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Because galactic cosmic rays also generate diffuse galactic radiation, it might appear difBcult to disentangle our 
signal from the background, which is dominant. Fortunately, it turns out that the energy spectrum of the photons has 
a well-defined peak, which is located approximately at the mass of the axion or ALP. This is expected to be of order 
/ieV, while diffuse galactic radiation tends to be much more energetic (GeV). The existence of such a peak can be 
understood as kinematics dictates that cosmic ray protons in the low energy end of their spectrum can only produce 
photons whose frequency is approximately the mass of the axion or ALP. Since both the cross section and cosmic-ray 
energy spectrum decrease with energy, photons whose frequency lies in the vicinity of the axion (or ALP) mass are 
most abundantly produced. 

The propagation of charged particles in a spatially homogeneous, but time-dependent, pseudoscalar background 
has also been investigated in |21| . in which the authors make the assumption that the time- varying background be 
treated as a constant in the Lagrangian (this is essentially the Lorentz-violating Chern-Simons term considered by 
Carroll et. al. in [22])- This additional term has the effect of modifying the dispersion relation of the photons. As 
a consequence, the process of photon emission is described at leading order by a single Feynman diagram with one 
vertex. According to |21j . photon emission by charged particles is then possible only if the photon is massive, and 
that the emission angle is small. In contrast, our calculation incorporates the time dependence of the axion field, 
which leads to three Feynman diagrams with two vertices each. Using this method, we find instead that emission is 
possible at all angles and for massless photons. 

This paper is structured as follows. In Section II, we review the physics behind the emission of photons by charged 
fermions propagating in a time-dependent pseudoscalar field. From this, we calculate the emission rate and energy 
spectrum of the photons. In Section III, we discuss the possibility of their detection with the Square Kilometre Array. 
In the absence of a detection, we put constraint on the value of fa- This is followed by conclusions in Section IV. 

II. THEORETICAL ANALYSIS OF PHOTON EMISSION BY PROTONS IN A PSEUDOSCALAR 

FIELD 

The Lagrangian that describes the dynamics of protons {ip) propagating in an axion (or ALP) field ((/)) is given by 

= -igapcjJtpJ^i^ + '^gajjFf.^F^"' + Cqed, (1) 

where Cqed is the usual QED Lagrangian that describes protons and photons. The couplings gapp = CapprUp/ fa 
and ga-yy — Ca-yyO:/ {2t: fa) are respectively the axion-proton and axion-photon coupling, where Capp and Ca-yy are 
dimensional-less model-dependent parameters, typically of order unity i23j. In this paper, we assume that Capp — 
Cay-y = 1. The parameter fa essentially measures the strength of the axion's coupling to the photon and proton. For 
the QCD axion, fa is known as the axion decay constant, and it is constrained to 10^ < fa < 10^^ GcV by particle 
and nuclear physics experiments, stellar evolution, and cosmology. This limits the QCD axion mass to [24 

« 6 X 10-6eV " . (2) 

Ja 

For ALPs, no such relation holds, and so both fa and nia are free parameters of the theory. 

The axion field is expected to oscillate, given by (/)(t) = 0o cos{mat) |24|- Because of its time variation and spatial 
homogeneity, it can be viewed as a source of energy, but not three-momentum. This is what makes the ordinarily 
kinematically forbidden process of photon emission by a charged particle possible. 
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A. Matrix elements 



There are three processes that contribute to the emission of photons by cosmic ray photons, and their corresponding 
Fcynman diagrams are shown in figures [T] [2] and [Sj The four-momenta of the incoming proton and the axion are 
denoted by and a'^, while those of the outgoing proton and emitted photon are and l'^ (i.e. p{p'^) + a{a^) — >■ 
p{q^) + 7(Z'^)). In the foUowing, we adopt the (H ) metric convention. 

The matrix element for the process ( "Process I" ; figure [T]) in which a cosmic ray proton emits a virtual photon, 
which then interacts with the axion field to emit a real photon, is 

= ^^^n{q)ru{p)e\P''{p~q)j0e;{l). 

The matrix element for the process ("Process 11"; figure [2]) in which a proton, subsequent to emitting a real photon, 
interacts with the axion field, is 

- ["(9)7'(2p^ - YlurXp)] ■ (3) 

The matrix element for the process ( "Process III" ; figure [s]) in which a proton interacts first with the axion 
background, then emits a real photon, is given by 

^M, = [u{q){2q^ + I^YDiMp)] ■ (4) 

To calculate the differential cross section, we need to first compute |A^2p, |A^3p, A^iA^2: MiMt,, ^^2-^3- 

Squaring Ali, averaging over initial proton spins, and summing over final photon polarizations, we have 

i^lXiP = -^^^[Aml{q-l){p-l)-A{p-q){p-l){q-l)+Arnl{l-af~2{l-af{p-q)]. (5) 

spins 

Meanwhile, squaring A^2 and yields 

2 E = 7^ [{P ■ l){q ■ I) + rnlia ■ I) - m^p ■ q) + m^] , (6) 

spins 

and 

2 E = 7^ [{P ■ l){q ■ I) + mlia ■ I) - m^p ■ q) + m^] . (7) 

spins 

The cross terms can likewise be straightforwardly computed: 

(^p)(p-9)2 

Ki^s il-q){p-qr 

lyM.Ml = -^^^[il.p)il.q) + {p.q){l.a)-ip-qf+mlip-q)]. (10) 

spins 



B. Differential cross section 



For simplicity, we evaluate the differential cross section in the rest frame of the axion field. In this case the four- 
momenta are given by p'' = {Ep, 0, 0,p), — {nia, 0), = (Eq, q sin 9, 0, q cos 9), = (w, Q), where Ep ~ ^p^ + m^, 

Eq — \ q^ + nip, and without loss of generality we align the z-axis with the direction of propagation of the initial 
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proton, and restrict the scattering to the x-z plane. 9 thus denotes the angle between the direction of the emitted 
photon and the z-axis. The photon's frequency uj is 



w? + 2maEp 

LU = . 

2Ep + 2ma - 2 |p I cose 

In this frame, the phase space for the final state particles is given by 



{2T:f2Eq (27r)32w 
which yields the differential cross section 



{2T:)H'^^\p + a-q~l), 



da 



iria + 2Ep 



dcos(e) iQTr{2Ep + 2ma - 2\p\ cos 0)2 /i^ 



(11) 



(12) 



C. Emission rate of the photons 



Our galaxy is teeming with cosmic rays, whose primary constituent is protons, to which we will restrict our attention 
in this paper (hence E = Ep below) . Our calculation of the photon emission rate is thus a conservative estimate, as 
other charged constituents (e.g. electrons) would also contribute to the process. As cosmic ray protons propagate 
in this background time-dependent axion field, they undergo photon emission via processes described in the previous 
section. 

To estimate the photon flux on Earth, we make the assumption that cosmic ray protons are homogeneous and 
isotropic within our galaxy. This is predicated on the observation that cosmic ray protons scatter off interstellar 
medium and traverse random trajectories within the Galaxy for an average of 10'' years in the galaxy j21j . For our 
calculation, we adopt the following energy spectrum for cosmic ray protons |26l I27j : 



dJ" / 3.06 \ f E \ 



dEdn \cm^ s sr GeV J \GeV J 



(13) 



which we assume to hold for i? > 50 GeV. Since the flux is known with less certainty for low-energy protons, we 
impose a cutoff in our calculation and disregard all protons with an energy below 50 GeV as a conservative measure. 
We will also neglect the contribution of extragalactic cosmic rays, which should be subdominant as compared to the 
galactic ones. 

Consider now a photon detector on Earth, with a field of view SQd- With our assumptions on the distribution of 
cosmic rays, the detector can pick up photons originating from cosmic rays filling a region from the Earth to the edge 
of the Galaxy, whose boundary is determined by 6ild- Let V denote this region, and dV = r^drSfld {r is the distance 
from dV to the detector on Earth) be a differential volume element in it. At each point in V, there are cosmic rays 
propagating in all possible directions. Locally, we define a spherical coordinate system (with the usual coordinates 9 
and 4>), where 9, as we defined earlier, denotes the angle between the velocity vector of the cosmic ray and the line 
connecting dV to the detector on Earth. The other variable (j) is the usual angle confined to the plane perpendicular 
to the direction of propagation of the cosmic ray. 

The number of photons emitted by the cosmic rays that fill up dV over an interval dt is given by 

dNj = UaUpdV daSvdt. (14) 

where Ua and Up are the number density of axions and protons, 5v the velocity of the cosmic rays, and da the 
differential cross section. From this, we can compute the the flux of photons per unit time at the detector, so we have 

'^'^^ ^J.dV, (15) 



'7 

dtdA \dE"" J dA 
where dA = r^dU — rodeos 9d4) is the differential area at the detector. Using the fact that 



dup ^ 1 dT 

dEdn ~ Sv dEdn ' ^ ' 
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the photon flux simphfies to 



dtdA 



Ua dFp da 

7^ dEdnd cose 



dVdEd COS 9. 



(17) 



Integrating over the volume that the detector can see, the proton energy, and cosO, and using (13 1 we obtain the 
photon flux at the detector. 



dN^ 



1 



dtdA 

which can be computed numericaUy. 



^dV ] I dE 



3.06 



i2 s sr GcV VGcVy 



E \ 



-2.70 



^ da{E,cost 
.1 dcos9 



-dcos t 



(18) 



D. Energy spectrum of photons 

Due to the uncertainties inherent in the energy spectrum of low-energy cosmic rays within the Galaxy, it is not 
possible to determine the precise spectral shape of the emitted photons. Nonetheless, the spectrum possesses a 
robust feature: it has a peak located at approximately the mass of the axion or ALP (more accurately, at lOc = 



(to^ + 2mamp) / (2mp + 2ma), obtained by setting Ep = nip in ( 11 )). 

The existence of such a peak can be understood as follows. From (11), we observe that for a given E, the emitted 
photon's frequency is conflned to the range cli_|- < lj < lo^, where 

^ — + (19) 

2E + 2maT'2x E'^ - ml 



The frequencies 0;+ and ui- correspond to scattering a.t 6 = and n respectively (see figure |4]). Hence, photons of 
frequency u can only be produced by cosmic ray protons with an energy larger than 

E = ( —iTL^ + 3m^uj — 2'm„uP' + \ l —Am^m'^^uj'^ + mf,uo'^ + 8m'},m„uj^ — 'im'^uj^ + 4m^ti;^^ . (20) 

2ml ~ AmaljJ \ y P o, a p a a j x i 

The frequency at which dE/duj vanishes is lj — uj^ (equivalently, E — mp). This implies that photons whose frequency 
is in the vicinity of Uc are produced by cosmic rays over the widest range of energy. More specifically, photons with a 
frequency near uj^. can be produced as long as the proton energy is above its rest mass. On the other hand, photons 
with a very low frequency near ma/ 2 can only be produced by enormously energetic protons which back-scatter at 
6 = TT, while those with a frequency oj ^ can only be produced by protons whose energy exceeds approximately 

^ojmp/{2ma)- Coupled to the fact that the cosmic ray flux and cross section decrease monotonically with E, we thus 

conclude that a peak is present in the energy spectrum of the photons. This has been verified numerically; see figure 
|6]for a plot of the spectrum, which features a salient peak, as expected. 



III. OBSERVATIONAL CONSEQUENCES 
A. Insufficient sensitivity to detect the QCD axion 



We numerically integrate ( 18 ) and find that, for = 10 ^ eV (corresponding to roughly the expected mass of the 
QCD axion), 

^ 10-2icm-2s-i. (21) 
dtdA ^ ' 

This estimate is quite conservative, as we only included energy £^ > 50 GeV in the integration, due to uncertainties in 
the spectrum of the protons. If we extend it to as low as, say, 10 GeV, we gain a boost in the fiux by approximately a 
factor of ten. From various measurements, we know that the flux actually increases as E is decreased to approximately 
1.4 GeV. For ria, we use Pa/ma, where pa = lO^^'' g/cm'^, the expected local halo density of the Galaxy [5S]. For the 
region V, we adopt the value 10^^ cm, which is approximately one tenth of the size of the stellar disk. The collecting 
area and field of view of the photon detector are taken to be that of the Square Kilometre Array: 10^° cm^ and 200 
deg^. Even with such a huge surface area, we expect only one photon every 10^^ s. Clearly, it is not yet possible to 
detect the QCD axion through this observation. 
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FIG. 4: Photon frequency u> versus the proton energy E, for an axion mass of 10~^ eV, according to ( |19[ ). To the right of the 
curve is the region in which photon generation by cosmic ray protons is aUowed. The turnaround point corresponds to u) = ijjc, 
which is very near the axion or ALP mass. The top (purple) and bottom (blue) curves correspond respectively to forward and 
backward scattering. 



dtr/dwx(dF/dE) 




FIG. 5: Photon energy spectrum dJ-/dE x da/duj versus photon energy lu for an axion mass of 10 ^ eV, up to a normahzation 
factor. For _E > 50 GeV, the proton spectrum is given by 3.06(£;/GeV)"^-^°. For nip < E < 50 GeV, it is chosen to be a linear 
function that extrapolates to zero a.t E = nip. This underestimates the abundance of low-energy cosmic ray protons, which 
however does capture the right picture that the flux vanishes at very low proton velocities. Note that in our actual calculation 
we make the even-more-conservative approximation of neglecting all cosmic rays whose E < 50 GeV. 



B. Constraining the parameter space of ALPs 

Although the detection of the QCD axion seems out of reach given existing technologies, interesting constraint can 
still be placed on the (ma, fa^) parameter space of ALPs, under the assumption that they constitute dark matter. 
Their number density is thus taken to be Pa/nia, where pa ^ 10^^'* g cm^'^. The sensitivity of the detector is 
assumed to be ten nanojanskys. Over a frequency range of 10^^ Hz, which is roughly what we need to see the peak in 
the spectrum, this translates to a minimum rate of 10"'^ photons/cm^s (for the Square Kilometre Array, about lO'' 
photons per second). Note that the constraint is only valid for an ALP whose mass lies between 2.9 x 10~^ eV and 
4.1 X 10^^ eV, since the Square Kilometre Array can only detect photons in this range with a field of view of 200 
deg'^. For higher photon energies (in the range 4.1 x 10~^ eV < rua < 1.2 x 10"'' eV that it is capable of detecting), 
the field of view diminishes substantially, resulting in a weaker exclusion limit. 

Our numerical result is shown in figure|6] The exclusion limit is slightly worse than that from laser experiments, but 
substantially weaker than that obtained from stellar evolution. It is foreseeable that improvements in radio telescopy 
in the future could help boost the sensitivity and thus improve the constraint. If so, our proposal would have the 
potential to complement existing axion experiments in the future, since it relies on very different physics from that 
in photon regeneration experiments (in fact, all other methods of axion search). 

Finally, we stress that the constraint shown in figure [6] is only approximate, given the uncertainties in the halo 
density, and the fiux spectrum and exact distribution of cosmic ray protons in the Galaxy. Throughout our calculation, 
we always try to be conservative in our estimates. For instance, we 1) disregard the contribution of all cosmic ray 
protons below 50 GeV, 2) consider only proton cosmic rays, and 3) only take into account cosmic rays within a 
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FIG. 6: Hypothetical exclusion limit (blue) on the mass ma and the decay constant fa for axion-like particles over mass range 
2.9 X 10"'^ eV <ma < 4.1 x 10"^ eV, in the case of a non-detection by the Square Kilometre Array. The other limits come 
from laser experiments, solar searches, and stellar evolution. The bottom line corresponds to the QCD axion. 

distance of 10^^ cm (one-tenth of the radius of the stellar disk) from Earth. 

IV. CONCLUSIONS 

In this paper, -we propose a ne-w observational probe of axions and ALPs, based on the detection of photons emitted 
by cosmic ray protons -which propagate in a time-dependent axion field -within the Galaxy. Since the axions couple 
very -weakly to ordinary matter, we find, unsurprisingly, that the photon emission rate is exceedingly lo-w, of order 
10~^^ cm~^ s^^. Even -with a detector -whose collecting area is as large as 10^" cm^ (e.g. the Square Kilometre 
Array), we expect only approximately one photon every 10^^ s. Thus detection of the QCD axion this -way is out of 
the question at the moment. 

Nonetheless, the same mechanism can be exploited to impose exclusion limits on the parameter space of ALPs. 
Since their decay constant fa can be smaller, ALPs can couple more strongly to ordinary matter, thereby increasing 
the rate of photon emission. Under the assumptions that dark matter is primarily ALPs, and that a detection rate of 
10^ photons per day is sufficient, -we numerically find that a non-detection of photons by the Square Kilometre Array 
translates to exclusion limits on fa (for the mass range 2.9 x 10~^ eV < uia < 4.1 x 10~^ eV) -which are slightly -worse 
than those obtained by current photon regeneration experiments. However, since the physics underlying our proposed 
observation is quite different from that used in any other axion/ ALP searches, it holds the promise for serving as an 
excellent complement to the various axion/ALP detection experiments in the future. 
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